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We present initial results from a wide field, multi-colour imaging project, designed 
to study galaxy evolution in X-ray selected clusters at intermediate (z ~ 0.25) and 
high redshifts (z ~ 0.5). In this paper we give blue galaxy fractions from eight X-ray 
selected clusters, drawn from a combined sample of three X-ray surveys. We find that 
all the clusters exhibit excess blue galaxy populations over the numbers observed in 
local systems, though a large scatter is present in the results. We find no significant 
correlation of blue fraction with redshift at z>0.2 although the large scatter could 
mask a positive trend. We also find no systematic trend of blue fraction with X-ray 
luminosity. We show that the blue fraction is a function of (a) radius within a cluster, 
(b) absolute magnitude and (c) the passbands used to measure the colour. We find 
that our blue fractions (ft,) from galaxy colours close to restframe (U — B)q, fb ~ 0.4, 
are systematically higher than those from restframe (B — V)o colours, fb ~ 0.2. We 
conclude this effect is real, may offer a partial explanation of the widely differing 
levels of blue fraction found in previous studies and may have implications for biases 
in optical samples selected in different bands. 

While the increasing blue fraction with radius can be interpreted as evidence of 
cluster infall of field galaxies, the exact physical processes which these galaxies undergo 
is unclear. We estimate that, in the cores of the more massive clusters, galaxies should 
be experiencing ram-pressure stripping of galactic gas by the intra-cluster medium. 
The fact that our low X-ray luminosity systems show a similar blue fraction as the 
high luminosity systems, as well as a significant blue fraction gradient with radius, 
implies other physical effects are also important. 

Key words: galaxies: clusters - galaxies: evolution - galaxies: fundamental param- 
eters - galaxies: elliptical and lenticular - galaxies: spiral 



1 INTRODUCTION 

During the past few decades, numerous studies of the op- 
tical properties of galaxy populations within clusters have 
allowed an intriguing insight into the processes that gov- 
ern the evolution of galaxies. Rich galaxy clusters provide a 
large sample of galaxies, at a similar redshift and of various 
masses and morphological type with which to probe these 
phenomena. One of the most influential early results, which 
motivated many successive studies was that of Butcher & 
Oemler (1978a). In this study of two high redshift clusters 



they found a large fraction of galaxies in the inner regions of 
the cluster with colours significantly bluer than the colour- 
magnitude relation (CMR) of the cluster ellipticals (Es). 
Butcher & Oemler (1984) extended this study to include 33 
clusters and found an increasing trend in this blue fraction 
(/s) with redshift, above z ~ 0.1. A number of subsequent 
photometric studies have generally confirmed the presence 
of these excess blue galaxies, though in varying quantities 
(e.g. Rakos & Schombert 1995; Smail et al. 1998; Margoniner 
& de Carvalho 2000; Kodama & Bower 2001). Subsequently 
spectroscopic investigations not only confirmed the cluster 
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membership of large numbers of these galaxies but also dis- 
covered that their spectra contained features indicative of 
star formation, often recently truncated (Dressier et al. 1985; 
Lavery & Henry 1986; Couch & Sharpies 1987; Dressier & 
Gunn 1992). 

During these studies it became increasingly clear that 
whilst the luminous elliptical galaxies present in cluster cores 
have undergone little but passive evolution since higher red- 
shift (z > 1, e.g. Ellis et al. 1997; Stanford et al. 1998), a ma- 
jor change has taken place in the other morphological types. 
Specifically the large SO population found in local clusters 
are found to be almost absent in high redshift (z ~ 0.5) clus- 
ters (Dressier et al. 1997), where evidence pointed towards 
the excess blue galaxies often having spiral morphologies. 
Confirmation of this came with morphological studies using 
HST. Dressier et al. (1994) and Couch et al. (1994) found 
that not only were the blue galaxy populations disk dom- 
inated, but they often exhibited signs of interaction. It is 
now generally thought that these photometric, spectroscopic 
and morphological evolutionary effects are related. The sug- 
gestion being that at higher redshifts, bluer field galaxies, 
generally spirals, that are drawn into the core regions of 
clusters have their star formation quenched by some mech- 
anism. They then proceed to evolve morphologies, possibly 
by some form of mass loss, into the large SO population seen 
in the core of rich clusters today. Perhaps the most surpris- 
ing aspect of these phenomena, though, is the rapid rate at 
which this evolution takes place. 

Clearly understanding the processes behind this evo- 
lution is vital to our understanding of how galaxy histo- 
ries have progressed. Whilst observational data is becoming 
ever more available, the mechanisms causing this evolution 
are less well understood. Several suggestions have been put 
forward to explain the effects of the interaction of galax- 
ies with their environment and each other. Certainly in the 
cores of clusters, where the intra-cluster medium (ICM) is 
densest, infalling galaxies, with smaller gravitational poten- 
tial wells, will suffer from ram-pressure stripping of their gas 
(e.g. Gunn & Gott 1972; Abadi et al. 1999). Other mecha- 
nisms which may affect smaller galaxies include tidal effects 
generated by the cluster potential (Byrd & Valtonen 1990) 
or galaxy harassment caused by close encounters with other 
galaxies (Lavery & Henry 1988; Moore et al. 1996). In ad- 
dition to these processes the direct interaction and merging 
of galaxies is also a possibility, especially in poorer clusters, 
with shallower gravitational potential wells (e.g. Toomre & 
Toomre 1972). 

Most studies till now have looked at massive clusters, 
nearly all of which have been optically selected. Several large 
photometric surveys, for instance, have used optically se- 
lected, such as Abell (1958), clusters (e.g. Butcher & Oemler 
1984; Rakos & Schombert 1995; Margoniner & de Carvalho 
2000). Differing blue fractions, however, were found in these 
investigations, with Rakos & Schombert (1995) claiming a 
blue fraction of up to 80 percent in high redshift clusters. 
Investigations based on clusters selected in other wavebands 
have only recently been undertaken. Allington-Smith et al. 
(1993) conducted a photometric study of galaxies around 
bright radio-galaxies up to z ~ 0.45, and found evidence 
for enhanced blue populations in most systems at higher 
redshift. The galaxy systems observed contained a number 
of poorer clusters and no relationship between blue galaxy 



fraction and redshift was found. They concluded that the 
increase in blue fraction in high redshift, rich clusters was 
due to processes inefficient in poorer systems. 

The other area which has only recently begun to be 
explored is the selection of clusters, for optical study, via 
the X-ray emission from their intra-cluster medium (ICM). 
Smail et al. (1998) used a sample of 10 intermediate redshift 
(z ~ 0.25), high X-ray luminosity, X-ray selected clusters 
in a multi-colour study of their galaxy populations. They 
found low blue fractions, with the bluer galaxies tending 
to avoid the cluster cores. The other main exploration of 
optical galaxy properties within X-ray selected clusters is 
that of the Canadian Network for Observational Cosmol- 
ogy (CNOC, Yee et al. 1996). They find excess blue galaxies 
in distant clusters (e.g. Morris et al. 1998), colour, emis- 
sion line, population and age gradients within the clusters 
(Balogh et al. 1997, Ellingson et al. 2000), and a fraction 
of luminous "K+A" cluster galaxies (with strong Balmer 
absorption lines but no emission lines) which is similar to 
that in the field at z^0.3 (Balogh et al. 1999). The com- 
bined N-body and semi-analytic models of Diaferio et al. 
(2001) match the CNOC1 cluster observations and support 
a gradual termination of star formation after cluster infall. 
Kodama & Bower (2001), using 7 clusters drawn from the 
CNOC sample also find an excess of blue galaxies. 

Whilst a general consensus as to the presence of this 
excess star-formation is evident, the level of it, and the evo- 
lutionary nature is less clear, with a major question being 
the steepness of the increase of blue galaxy population with 
redshift (e.g. Margoniner & de Carvalho 2000) . Blue fraction 
may, or may not, be a function of redshift and X-ray lumi- 
nosity, which may cause significant biases in studies limited 
to a small sample of clusters, within a limited range of red- 
shifts or samples compiled for comparison which have vary- 
ing X-ray luminosities (e.g. Andreon & Ettori 1999). The 
works of Lea & Henry (1988), Smail et al. (1998) and An- 
dreon & Ettori (1999) hint at there being no clear relation 
between blue fraction and X-ray luminosity. Additionally, 
the changing environmental nature of the clusters may also 
have a large effect, for instance excess sub-structure may 
bias towards higher blue fractions (Metevier et al. 2000). It 
can, therefore, be seen that the large spread in observed re- 
sults highlights the need for a method of cluster selection 
which avoids the biases inherent in optical selection, whilst 
still providing large catalogues of galaxy clusters with multi- 
farious properties to allow exploration of the redshift, cluster 
properties and environment parameter space. 

This paper, the first in a series, presents details of anal- 
ysis and initial results from multi-colour, wide-field imaging 
of an X-ray selected sample of galaxy clusters. The clusters 
span a large range in X-ray luminosity (a factor of 100 in 
the 8 clusters presented here) whilst being constrained to 
two redshift bands designed to minimise the effects of k- 
correction. This project, on completion, should allow an in- 
sight into the evolution of galaxies within a variety of cluster 
environments, selected solely according to their X-ray prop- 
erties. The organisation of the paper is as follows. In the next 
section we discuss further the subject of the selection of the 
sample of clusters presented here. We then detail the obser- 
vations used in this paper, before discussing the reduction 
and analysis techniques used. We then present and discuss 
results from the photometry, including colour-magnitude di- 
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agrams and blue fractions. Our conclusions are drawn in the 
final section. In order to allow comparison to other works we 
consider an Ho = 50kms _1 Mpc - and qo = O.f cosmology 
throughout this paper. 



2 SAMPLE SELECTION 

As discussed in the previous section, the need for an unbi- 
ased cluster sample with which to investigate galaxy evolu- 
tion is of paramount importance. This is certainly the case 
when trying to construct a picture of the various mecha- 
nisms that affect the changing morphologies and star for- 
mation histories of cluster galaxies. Whilst the studies of 
optically rich clusters have given evidence for rapid evolution 
in galaxy colours, this is by no means a universal result. In- 
deed even within the Butcher & Oemler (1984) analysis, the 
scatter on the blue fraction versus redshift plot is large. For 
instance CL0016+16, A2397 and A2645 (at redshifts 0.541, 
0.224 and 0.246 respectively), to name but a few, all demon- 
strate little, if any, evidence for enhanced blue fraction. This 
may be indicative of the fact that optically selected high 
redshift clusters, which traditionally were often selected in 
a blue rest frame (e.g. Abell 1958; Gunn et al. 1986), may 
be an unrepresentative sample and biased towards high blue 
fractions. Indeed Smail et al. (1998), using an X-ray selected 
sample of X-ray luminous massive clusters found evidence of 
a spread in blue fractions, though with a low median value 
(fb = 0.04 for concentrated clusters, for which Butcher & 
Oemler (1984) found fb = 0.09 for a similar sample), and 
hence generally low star formation rates. The goal of this 
project, therefore, is to construct a sample of galaxy clus- 
ters where the selection of the clusters is independent of the 
galaxy properties, with which to investigate cluster galaxy 
evolution. Additionally we wish to explore a wide range of 
cluster masses (as indicated by their X-ray luminosities) and 
be able to make comparisons between similar mass clusters 
at different redshifts. 

The clusters used in this study are drawn from the 
catalogues of three major cluster surveys. The Einstein 
Medium Sensitivity Survey (EMSS, Gioia & Luppino 1994) 
has a higher flux limit than the two other surveys used and 
thus primarily contributes brighter X-ray clusters. In addi- 
tion we use two other serendipitous cluster catalogues based 
on archival, ROSAT PSPC data. These are the Southern- 
Serendipitous High-redshift Archival ROSAT Cluster sur- 
vey (S-SHARC, Burke et al. 1997; Collins et al. 1997) and 
the Wide Angle ROSAT Pointed Survey (WARPS, Scharf 
et al. 1997; Jones et al. 1998) . All of these surveys use the X- 
ray emission from diffuse gas, trapped in the gravitational 
potential well of a cluster, as tracers for the mass of the 
cluster. Optical imaging and spectroscopic follow-up then 
confirms the presence and redshift of the clusters from a 
few of their brightest galaxies. This method has the consid- 
erable advantage over optical surveys of nearly eliminating 
biasing caused by projection effects and erroneous identi- 
fications (Lucey 1983; Frenk et al. 1990; Struble & Rood 
1991). Whilst the three X-ray surveys used as a basis for 
this study vary in their flux limits and detection techniques, 
they are all based on the detection of an X-ray flux in a 
similar energy range, and they all corrected from measured 
luminosities to total luminosities in a similar manner, as- 



suming a King surface brightness profile. Mean differences 
in total flux estimates between X-ray surveys are typically 
f»20 percent (e.g. Jones et al. 1998, Vikhlinin et al. 1998). 

We attempted to restrain the redshift ranges of our 
cluster sample to minimise the effects of k-correction be- 
tween bands. To this end we selected two redshift bands for 
our clusters, in which redshift corrected standard Johnson- 
Kron-Cousins filters (BVR in the low redshift sample, VRI 
in the high) were mapped onto rest frame UBV filters, as 
much as possible. The limited number of X-ray selected clus- 
ters in any given RA range, however, prevents a very tight 
constraint in redshift range. For the purposes of this initial 
study, eight clusters were included in our sample for obser- 
vation. These were chosen to maximise the spread in X-ray 
luminosity within the constraints placed by the RA range 
of this first observing run. Details of these clusters can be 
found in Table |l[ Here we give cluster IDs along with spec- 
troscopically confirmed cluster redshifts. All the clusters in 
this sample have a reasonably relaxed X-ray morphology as 
detected in their ROSAT PSPC data. They should, thus, 
be good examples of reasonably relaxed systems. The blue 
fraction of one of our sample, MS1621.5, has been studied 
previously (Morris et al. 1998; Kodama & Bower 2001) and 
thus will allow comparison to this investigation, providing 
a good consistency check. The sample presented here is by 
no means complete. Future papers will detail observations 
expanding the sample to fill out the redshift and X-ray lu- 
minosity parameter space. 



3 OBSERVATIONS 

The eight clusters detailed in Table [j] were observed on June 
17-20th 1999 using the Wide Field Camera (WFC, Ives et al. 
1996) on the 2.5 metre Isaac Newton Telescope (INT) on La 
Palma. The WFC consists of 4 thinned EEV CCDs each 
having 0.33 arcsec/pixel resolution and giving a combined 
spatial coverage of ~ 0.29 deg 2 (each 22.8 by 11.4 arcmin) 
at telescope prime focus. The 4 CCDs exhibit varying char- 
acteristics with, in particular, the gain of the central CCD 
being lower than the others. In addition there is some non- 
linearity apparent in all CCDs, with the most non-linear 
varying by around 6 percent over the entire dynamic range. 
The large sky coverage of the WFC allows us, generally, to 
image the entirety of the target cluster in the central CCD. 
In the cosmology used in this paper, Ho = 50kms _1 Mpc -1 
and qo = 0.1, the 5.7 arcmin half width of our central CCD 
corresponds to a little over 1.7 Mpc, for our nearest cluster, 
and around 2.9 Mpc for our most distant. In addition the 
WFC also provides a highly accurate statistical estimation 
of the local background of the cluster, without the need for 
separate calibration fields. This will hopefully reduce the af- 
fect of large-scale structural variations on our investigation. 
We observed each of the eight clusters in three filters, BVR 
and VRI for the low and high redshift clusters respectively. 
Details of the exposure times obtained are given in Table |l| 
In each case we aimed to achieve good signal-to-noise down 
to at least M* + 3 (Mv ~ —19). In addition we give the clus- 
ter X-ray luminosity from the ROSAT PSPC (0.5-2 keV). 
Most observations were carried out near zenith and seeing 
conditions were found to be generally good, ranging from 
1.52 to 0.79 arcsec FWHM. 
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Cluster 


Rcdshift 


L x (crgs 1 ) 


Filters 


Exposure (ks) 


Completeness 


RXJ1633.6+5714 


0.239 


0.12 x 10 44 


B,V,R 


3.00, 1.20, 1.20 


24.3, 24.2, 23.5 


MS1455. 0+2232 


0.259 


11.51 x 10 44 


B, V, R 


2.24, 2.40, 2.70 


24.0, 23.9, 22.5 


RXJ1418.5+2510 


0.294 


3.66 x 10 44 


B,V,R 


3.60, 1.20, 1.20 


24.2, 23.9, 22.4 


RXJ1606.7+2329 


0.310 


0.85 x 10 44 


B,V,R 


3.60, 1.20, 1.20 


24.2, 24.1, 22.5 


MS1621. 5+2640 


0.426 


4.71 x 10 44 


V,R,I 


8.40, 5.40, 1.80 


24.0, 23.9, 22.9 


RXJ2106.8-0510 


0.449 


2.79 x 10 44 


V,R,I 


6.90, 4.20, 1.80 


24.2, 24.1, 22.9 


RXJ2146.0+0423 


0.532 


4.27 x 10 44 


V,R,I 


6.00, 2.40, 1.80 


24.4, 24.2, 22.8 


MS2053. 7-0449 


0.583 


5.49 x 10 44 


V,R,I 


6.00, 3.90, 1.80 


24.0, 23.8, 22.8 



Table 1. Details of the 8 X-ray selected clusters. Redshifts are spectroscopically confirmed. X-ray luminosities from 
the ROSAT PSPC are given in the 0.5-2 keV band. Listed are the filters used, along with the respective exposure 
times. Finally apparent magnitude completeness limits are given, estimated from N(m) plots. 



4 DATA REDUCTION AND ANALYSIS 
4.1 Reduction and calibration 

Initial reduction was carried out in the standard way. Af- 
ter bias subtraction and non-linearity corrections (obtained 
from the Wide Field Survey (WFS, McMahon et al. 1999)) 
were applied, the images were flat -fielded. Large scale gradi- 
ents were observed in the B, V and R filters after application 
of twilight flats, so instead it was decided to use data flats 
median combined from the dithered science images, with 
masks designed to eliminate bright galaxies and stars. Due, 
however, to the low number of target clusters observed dur- 
ing each night, allied to lengthy exposures (chosen because 
of long readout times), the number of images available for 
median combination each night was insufficient to generate 
statistically acceptable flat-fields. It was therefore necessary 
to combine data from more than one night. We used median 
combined data flats combined from two nights, in the V 
and R bands. In the case of the B band exposures, data was 
taken on only three nights, and all were required to provide 
an adequate flat. In the cases of the I and R band expo- 
sures, night sky fringing was apparent. In the R band this 
was found to be virtually eliminated by flat-fielding using 
data flats. The level of fringing in the I band, however, was 
significantly higher (of the order of three percent deviation 
from sky level) . The reduction of the I band thus followed a 
slightly different pattern. Here we used the twilight flat-field 
which produced highly fringed, though generally flat images. 
We proceeded to produce a master fringe frame from all of 
our I band data and this was used to de-fringe our data. We 
used an iterative process to scale the master fringe frame to 
the approximate level of the fringes; subtracting this left a 
residual fringe level in each image, which could be used to 
adjust the subsequent level of the scaling. After final de- 
fringing, our images typically showed residual fringing at 
IS 0.3 percent of sky level. The above steps were carried out 
independently on each CCD. Figure |l| contains an example 
V band image of the central region of RX J 14 18. 5. Over- 
layed upon this image are X-ray contours from the ROSAT 
PSPC. 

At least twelve sets of Landolt (1992) standard star 
fields (Landolt 107,110,113) were observed in various filters, 
each night, to allow accurate photometric zero-points to be 
obtained for each CCD. Our first night was found to be 
non-photometric, though in all cases clusters observed dur- 
ing this night had short exposure calibration images taken 
on one of the following nights. Nights 3 and 4 of the ob- 



servations appeared to be totally photometric and the at- 
mospheric absorption on these nights was found to be con- 
sistent with the standard values of the La Palma site, and 
thus we adopted these. It was found that an excess amount 
of atmospheric absorption was present on the second night 
(R « 0.30 mag at zenith). This appeared to be present and 
contributing similar absorption in all bands and did not ap- 
pear to vary throughout the night. We attribute this roughly 
grey absorption to atmospheric dust. It was also observed 
(with a slightly larger magnitude (r 1 ~ 0.38)) in the r' band 
standard star measurements at the Carlsberg Meridian tele- 
scope, on La Palma, during the same period. Zero points 
and colour equations were devised for each CCD, in each 
filter on each night, in order to place the photometry on the 
standard Johnson-Kron-Cousins system. Although the final 
solutions to these equations had rms deviations of less than 
0.005 mag, the maximum final absolute photometric error, 
including systematic effects, was estimated to be < 0.03 
mag. This estimation was based on standard star deviations 
from our best fit solutions, as well as from external checks 
which are discussed in the next section. In general the so- 
lutions were consistent with those obtained by the Wide 
Field Survey. In addition zero-points were adjusted for the 
affects of galactic absorption, for each cluster pointing, us- 
ing absorption maps from the IRAS satellite (Schlegel et al. 
1998). 



4.2 Constructing photometric catalogues 

Galaxy photometry catalogues were constructed using the 
SExtractor package (Bertin & Arnouts 1996). This allows 
detection of objects in one image, with additional photo- 
metric analysis from another. Our images were aligned and 
trimmed to provide the same sky coverage for all filters for 
each CCD. Detection was carried out in the reddest filter 
image for each cluster (i.e. I band for the higher redshift 
clusters, R band for the lower) so as not to exclude the 
redder cluster galaxies. Detection significances of more than 
1.5 a per pixel were required consisting of at least 9 con- 
tiguous pixels. We used the observed colour of each galaxy, 
with the appropriate colour equation, to give accurate mag- 
nitudes. We used matching aperture magnitudes for colour 
estimation, with the aperture size, a 6 pixel radius equiva- 
lent to 2.0 arcsec, larger than the worst seeing in any of our 
images. Kron-type magnitudes were measured, using the de- 
fault SExtractor settings, to provide the best estimate of the 
total magnitudes of the galaxies. 
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Figure 2. A composite number-magnitude diagram for the combined backgrounds of the clusters in the Ft band. Different symbols 
represent the 8 different clusters. Error bars are included but are generally smaller than the point size. The lines overlayed are from 
the works of Jones et al. (1991); Metcalfe et al. (1991, 1995); Steidel & Hamilton (1993) and Smail et al. (1995). The turnover due to 
incompleteness generally occurs at brighter magnitudes for the lower redshift clusters. 



After elimination of objects with unreliable photome- 
try due to, for instance, proximity to CCD edge, star/galaxy 
separation was achieved using the SExtractor CLASS-STAR 
parameter estimated from the two reddest bands for each 
cluster. At fainter magnitudes all star/galaxy classification 
begins to fail, thus we limited our separation to brighter than 
an apparent magnitude of I = 22, R = 22.5 and V = 23. 
A number of fainter stars will thus be retained, however at 
these fainter magnitudes galaxy counts dominate (e.g. Met- 
calfe et al. 1991) and the stars will remain in both cluster 
and field images and so should be naturally accounted for. 
In the case of SExtractor detection in our / band images, a 
few false detections were noticed by visual inspection, these 
coincided with residual fringes in the images. The photom- 
etry extracted from the other two filters, however, showed 
no sign of these detections and so provided a useful robust- 
ness check. All galaxies not detected in at least two of the 
bands were excluded from our analysis. This cut may also 
exclude objects of extreme colour towards our magnitude 
limit, these however generally fall well below the magnitude 
cuts used in the analysis presented here. 

Number-magnitude plots of galaxies were generated al- 
lowing us to make conservative estimates of the apparent 
magnitude completeness level of each observation. An ex- 



ample of one of these plots is shown in Figure |2|. Here we 
plot the number-magnitude relations for the R band obser- 
vations of the combined background fields of all our clus- 
ters, in 0.5 magnitude bins. Overlayed on this plot are best 
fit lines from previous large scale galaxy number counts in- 
vestigations (Jones et al. 1991; Metcalfe at al. 1991, 1995; 
Steidel & Hamilton 1993; Smail et al. 1995). It can be seen 
that, whilst all of our backgrounds are in reasonable agree- 
ment with the previous studies, a significant spread in nor- 
malisation is found (up to a factor of around 1.5 difference). 
This is important as it emphasizes the advantage of hav- 
ing large background areas surrounding the cluster, which 
can then be used to take account of variations in the field 
galaxy density used for the statistical analysis. Complete- 
ness levels estimated from the number-magnitude plots are 
also given in Table [j]. The variations between exposure times 
and completeness levels are due to the varying observa- 
tional conditions. In most cases the completeness levels cor- 
respond to reaching our original target depth of My = — 19 
(one magnitude deeper than that used in Butcher & Oemler 
(1984)) in our cosmology. We miss My = —19 by about 0.5 
magnitudes for our most distant cluster, MS2053. 7-0449. In 
addition to the calibration provided by our standard star 
measurements we use two other independent diagnostics to 
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Figure 1. INT WFC V band image of RXJ1418. 5+2510 at z = 
0.294. Overlayed on the image are contours of the X-ray emission 
detected from this L x = 3.66 X lO^ergs -1 (0.5-2 keV) cluster 
from the ROSAT PSPC. Contours are equally spaced in flux level 
and were derived from an image smoothed to PSPC resolution. 
Only a small fraction (~ 2 percent) of the total WFC image is 
shown. 



confirm our photometry. The field of RXJ1418. 5+2510 has 
Johnson B-band magnitudes of several stellar sources pub- 
lished by Tyson et al. (1998), which agreed to B > 21.5 
with our measured magnitudes, to within 0.03 mag rms. The 
other checks are of a more statistical nature. We generated 
galaxy and stellar count and colour histograms for several 
of the background fields of our clusters. These were then 
compared to those published by Metcalfe et al. (1991) and 
again gave good agreement, given statistical fluctuations. 
For instance our average (B — R) CCD colours of field galax- 
ies (when converted to the filter colours of Metcalfe et al. 
(1991)) for our four low redshift cluster fields were 1.79, 
1.75 and 1.79 (in the magnitude ranges 19.5 < Rccd < 20, 
20 < Rccd < 20.5 and 20.5 < Rccd < 21 respectively) which 
compare well with the corresponding values of 1.70, 1.93 
and 1.78 (from Metcalfe et al. (1991)), given the sizeable 
statistical errors and large scale structure variations. 

4.3 Galaxy density profiles 

The variation in cluster mass built into our sample means 
that each cluster has a different physical size. This is impor- 
tant to take into consideration when comparing the popula- 
tions of the various clusters. Matched physical sizes will not 
provide an accurate comparison, as clusters cores are known 
to contain galaxy populations with different characteristics 
to those in more outlying regions (e.g. Dressier et al. 1997; 
Morris et al. 1998). A better analysis technique would be to 
use the galaxy population as an indicator of the cluster size. 

Butcher & Oemler (1978b) detail a method of estimat- 
ing matched radial distances in clusters from their galaxy 
density profiles. They define a concentration parameter 



C = log(R 6 o/R 2 



(1) 



where R n represents the radius containing n percent of 
the total cluster population. For each of our clusters, radial 
galaxy density profiles were generated, allowing estimates 
of the cluster concentration parameter and the radius of the 
cluster which contained 30 percent of the cluster galaxies 
(R30, as used for blue fraction calculations in Butcher & 
Oemler 1984). 

It was clear from the statistically background sub- 
tracted radial profiles where the cluster over-density reached 
the background levels, and this allowed estimates of the 
characteristic radii from the integrated population within 
this boundary. The large background fields available for each 
cluster served to reduce the error in the background level es- 
timates. The derived values of R30 were not sensitive to the 
details of this procedure, and for all our clusters R30 lay 
well within the cluster CCD. Our R30 and concentration 
parameter values for each cluster are detailed in Table ^[ It 
is reassuring to note that all our clusters have concentration 
parameters above 0.4, which would be expected from the re- 
laxed nature of their X-ray contours. It is these concentrated 
clusters which were used by Butcher & Oemler (1984) to fit 
a trend line which defined the Butcher-Oemler effect. 



4.4 Bright galaxy area correction 

The main complication of generating accurate radial pro- 
files is the estimation of the number of fainter galaxies lost 
behind brighter and larger galaxies in the foreground. To 
solve this problem we calculated the total area lost behind 
brighter galaxies, in each radial bin used for this calcula- 
tion. We chose all galaxies above 300 pixels in size and es- 
timated the area they occupied, within the radial bin area, 
at a brightness level above their detection threshold where 
blended galaxies would not have been detected. The area 
used in the calculations was then corrected for this lost 
area. On average the correction was of order one percent 
for the field frames and around two percent for the cluster 
fields, which agrees well with previous analyses (e.g. Tren- 
tham 1997). 



5 RESULTS 

5.1 Colour— magnitude diagrams 

Cluster plus field colour-magnitude (CM) diagrams were 
generated, along with field CM diagrams from the three 
non-cluster CCDs. This allowed a statistical background 
subtraction to be carried out, for the purposes of illustration 
of cluster against field galaxies. This was achieved by assign- 
ing a membership probability to each galaxy, determined 
from the expected number of galaxies from the background 
fields. We then used a Monte-Carlo technique, analogous to 
that used in Kodama & Bower (2001), to assign probabilities 
of cluster membership to each galaxy. The large background 
area meant that this could be achieved in small boxes of size 
0.1 mag in colour and 0.5 mag in apparent magnitude. 

In Figure ^| we present five statistically background sub- 
tracted CM diagrams, in which we show colours correspond- 
ing approximately to restframe B — V versus V. We plot the 
apparent 7? magnitude against (V — R) colour for MS1455.0, 
an X-ray luminous cluster and similarly we show (V — R) 
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Figure 3. Colour— magnitude diagrams for five X-ray selected clusters. Here we plot the colours which correspond approximately to 
(B - V)o versus V at rest: observed V - R against R for MS1455. 0+2232 (z = 0.259, L x = 11.51 X 10 44 ergs" 1 ), RXJ1633.6+5714 
(z = 0.239, L x = 0.12 X 10 44 ergs" 1 ) and RXJ1606.7+2329 (z = 0.310, L x = 0.854 X 10 44 ), and R - I against I for RXJ2106. 8-0510 
(z = 0.449, L x = 2.9 X 10 44 ergs" 1 ) and MS2053. 7-0449 (z = 0.583, L x = 5.5 X 10 44 ergs" 1 ). Also shown is a composite plot of the CM 
diagrams of the background fields of RXJ1606. Only galaxies within R30 are shown. Our best fitting CMRs arc shown as the solid lines. 
Blue galaxy cuts used in determining the blue fraction are shown as the dashed lines below the CMR line. Vertical dot-dashed lines 
represent a rest frame My = —20 magnitude cut, whilst the triple dotted-dashed line diverging from this illustrates the effect of the 
different level of k-corrections for varying galaxy colours. The small boxes at the bottom of each CM diagram represent average colour 
and magnitude photometry error boxes, at the corresponding magnitudes. Finally the plots have sloped dotted lines which illustrates 
if and where our 100 percent completeness level impinges upon the CM diagram. 
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Figure 4. Colour— magnitude diagrams for the three remaining X-ray selected clusters. Here we plot the colours which correspond 
approximately to (U - B) versus V at rest: observed V - R against / for MS1621. 5+2640 (z = 0.426, L x = 4.71 X 10 44 ergs -1 ) and 
RXJ2146.0+0423 (z = 0.532, L x = 4.27 X 10 44 erg s -1 ), and B — V against R for RXJ1418. 5+2510 (z = 0.294, L x = 3.66 X 10 44 erg s -1 ). 
Also shown is a composite plot of the CM diagrams of the background fields surrounding MS1621.5. Again in all these cases galaxies 
within R30 are shown. The lines drawn on the plots represent the same as in the previous figure. 



against 7? for RXJ1633.6 and for RXJ1606.7, two of our 
low L x clusters. We also plot apparent I magnitude against 
(R - I) colour for RXJ2106.8 and MS2053.7, two high red- 
shift clusters. Crosses represent galaxies statistically flagged 
as cluster members and dots represent the statistical field 
population. Additionally the corresponding composite field 
CM plot for the background fields of RXJ1606.7 is also plot- 
ted. 

In Figure ^ we show CM plots of the three remaining 
clusters. Here though we plot colours corresponding approx- 
imately to restframe U — B versus V. Additionally the corre- 
sponding composite field CM plot for the background fields 
of MS1621.5 is shown opposite it. At the bottom of each of 
these plots are boxes representing an average error on our 
photometry. The errors on the colours were taken from the 
error on the aperture photometry and added in quadrature. 
As can be easily seen, only in our faintest galaxies does the 



error on our colours begin to reach a level where it would 
effect our results. 

Also shown in Figure [| are our completeness limits, 
if they impinge upon the CM graph. In these cases, only 
RXJ2106.8 is affected. The sloped, dotted line in the top- 
right hand corner represents our conservative estimate of the 
100 percent completeness level. We do not expect this limit 
to seriously effect our blue fraction results, at the magni- 
tude levels of ~ Mv = —20. Additional evidence for this 
comes from the corresponding CM diagram of MS1621.5, 
from Kodama & Bower (2001), which does not show any 
excess red population near our completeness limit. In the 
cases of the graphs displayed in Figure ^, the complete- 
ness limit is harder to display, as it is a function of all 
three filter magnitudes. We make an attempt to constrain 
these by using colour-colour plots of our galaxies, which 
will be presented in more detail in a forthcoming paper. We 
find that for our low redshift sample we are 100 percent 
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complete in B — V at Mv = —20, and estimate that this 
drops to around 95 percent completeness at Mv = —19. 
For our higher redshift sample the results are poorer. Only 
in MS1621.5 are we fully complete to M v = -19 . In the 
cases of RXJ2106.8 and RXJ2146.0 we are around 95 per- 
cent complete at Mv = —20. We are significantly incomplete 
to M v = -20 for our furthest cluster, MS2053.7, and will 
thus not give results for the V — R blue fraction of this 
cluster. The effect of the incompleteness, in the cases where 
it may impinge on our CM diagrams, would be to increase 
the value of the blue fractions marginally. Our estimates of 
missing galaxies are confirmed by comparison with the ob- 
servations of Metcalfe et al. (1991) who present histogram 
of galaxy numbers versus colour. Our field samples are in 
excellent agreement with these histograms. 

5.2 Colour— magnitude relation 

The CMR of cluster E/SOs was fitted using a bi-weight 
estimator (Beers et al. 1990). In our two poorest clusters, 
(RXJ1633.6 and RXJ1606.7) we had too few cluster galax- 
ies to fit the CMR. Here we took advantage of the similarity 
in redshift and X-ray luminosity of these 2 systems. We com- 
bined the CM diagrams of these two systems, correcting for 
the difference in E/SO colour due to redshift (Fukugita et al. 
1995) and fit for the CMR from the combined CM diagram. 
In most cases these bi-weight fits provided what appeared 
to be a good estimate of the slope and normalisation of the 
CMR. In certain cases, generally those with poor statistics 
caused by lack of galaxies, the CMR fit was not considered 
to be an accurate reflection of the actual CMR. In these 
cases we took a 0.3 magnitude cut either side of the ex- 
pected E/S0 colour at R or J = 20, and refit the CMR with 
a weighted least squares fitting algorithm. In general this 
served to alter the normalisation of the CMR without alter- 
ing the slope significantly. The CMRs of the eight clusters 
all have a similar slope and are consistent with a universal 
CMR, given the errors on the fits. Detailed analysis of the 
cluster CMRs will be presented in a future paper. 

Finally, to avoid the problems inherent in k-correction 
of galaxy magnitudes to rest frame, we use the approach of 
Kodama & Bower (2001) and convert the Butcher & Oemler 
(1984) blue galaxy cut to the observed cluster redshift. This 
was done by taking the A(B — V)o ~ —0.2 colour cut from 
the CMR to approximately represent an Sb/Sbc galaxy and 
estimating, via interpolation between the different redshift 
and morphological values, the corresponding colour differ- 
ence at the cluster redshift using observed colours (Fukugita 
et al. 1995). Overlayed on the CM diagrams in Figure^ and 
Figure [l| are the fitted CM relations along with the result- 
ing blue galaxy colour cut. The graphs display a magnitude 
cut equivalent to M v = —20, as used in Butcher & Oemler 
(1984). This was calculated assuming average E/SO colours 
and in a cosmology with Ho = 50kms _1 Mpc _1 and qo = 0.1, 
and is shown as the vertical dashed-dotted line. This, how- 
ever, will tend to include or exclude bluer galaxies, as their k- 
corrections are colour dependent. We therefore use observed 
galaxy colours (Fukugita et al. 1995) to modify this magni- 
tude cut to account for this. These modifications are shown 
as the triple-dotted-dashed lines. Where the M v = — 20 line 
and the colour corrected magnitude cut line overlap almost 
completely (e.g in the CM plot of RXJ1418.5 or RXJ2106.8) 



we find that the observed filter maps extremely well onto the 
restframe filter for all galaxy colours. 



5.3 Blue fractions 

Whilst the statistical subtraction of background galaxies 
from the cluster plus field CM diagram gives a good indica- 
tion of the expected distribution of cluster galaxies within 
these CM plots, the accuracy is limited by the size of colour 
and magnitude bins used. This then may bias any calcula- 
tion involving the galaxies selected within these boxes. A 
better way to utilise the large background areas surround- 
ing our clusters would be to directly calculate the expected 
number of background blue versus red galaxies, and then 
adjust the cluster plus field results accordingly. This avoids 
any problems of colour-magnitude bins bisecting the CMR 
or blue cut line. This then was the method used for our anal- 
ysis. We used the largest background area available, at dis- 
tances always greater than 7.5 arcmin from the cluster cen- 
tre, to calculate the most statistically accurate background 
estimate possible. Again the area lost to brighter galaxies 
was taken into account during the scaling to the cluster area 
(see Section 4.4 ) . Simple subtraction of statistical field blue 
number from actual cluster plus field blue number, and back- 
ground total from cluster plus field total allowed an estimate 
of the blue fraction to be found. We assumed that the error 
on each of these terms could be added in quadrature. Un- 
surprisingly the dominating error on the final results comes 
from the relatively small excess of blue cluster galaxies. 

Table ^| shows blue fractions (fb), calculated at R30 and 
with a magnitude cut of M v = — 20 for our eight clusters. 
It should be noted that the V — R and R — I blue fractions 
for our higher redshift clusters, and the B — V and V — Rib 
values for our lower redshift clusters, roughly correspond to 
U — B and B — V respectively at rest. In addition to the blue 
fractions derived from the method detailed above, we also 
compare our estimates to those obtained from the method 
of Kodama & Bower (2001). Namely using the colour mag- 
nitude diagrams with statistical background subtraction in 
colour and magnitude bins. As can be seen from these fig- 
ures given in Table ^, we find a general agreement in results 
between the two methods, within errors. The other compar- 
ison available to us is, in the case of MS1621. 5+2640, with 
published values of fb. Whilst our results are not directly 
comparable to those of Morris et al. (1998) and Kodama & 
Bower (2001), due to the different filters used in the observa- 
tions, we find a generally good agreement. Firstly our value 
of R30, 2.24 arcmin, lies close to that of 2.18 arcmin as given 
in Kodama & Bower (2001). Secondly our blue fractions in 
both (V - R) and (R - J), 0.13 ± 0.04 and 0.22 ± 0.05 re- 
spectively, agree with the values of fb = 0.16 ± 0.04 and 
fb ~ 0.2 from the studies of Kodama & Bower (2001) and 
Morris et al. (1998) respectively. 



5.4 Individual cluster anomalies 

Here we briefly discuss issues arising from varying aspects 
of the differing properties of our clusters, as well as note any 
anomalies in our analysis procedure for individual clusters. 

RXJ2146. 0+0423: The CM diagram for this cluster 
(see Figure m shows some evidence of contamination from 
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Table 2. Details of the 8 X-ray selected clusters. X-ray luminosities from the ROSAT PSPC are given in the 0.5-2 
keV band. Concentration parameters (C) and R30 values are derived from background subtracted galaxy radial 
profiles. Blue fractions are given, for M v < —20, in B — V and V — R colours for the low rcdshift sample, and 
V — R and R — I for the higher rcdshift clusters. Additionally, for comparison, in the final two columns, we give 
blue fractions derived from CM diagrams field subtracted in colour and magnitude boxes (see text). 



another galaxy cluster system. This is evident in what may 
be a second CMR at around the V — 0.6 level. No X-ray 
contamination of this cluster is detected, hence the inclu- 
sion of it within our original sample. Additionally no obvi- 
ous galaxy excess is visible outside the target cluster. This 
is confirmed by the survey of Gunn et al. (1986) who detect 
our target cluster, but find no others within the area of our 
central CCD. We note, therefore, that the estimation of the 
blue fraction from this cluster is probably incorrectly high. 
We also believe this highlights the possibilities of serendip- 
itous group detection in cluster CM studies. In addition we 
note that similar effects are observed in other studies. For in- 
stance the CM diagram of MS1512.4+3647 within the study 
by Kodama & Bower (2001) displays a similar effect, which 
the authors attribute to foreground group contamination. 

MS1455. 0+2232: This cluster has the largest calcu- 
lated value of R30 in our sample (3.75 arcmin) and whilst 
this radius lies comfortably within the central CCD, the ra- 
dial density profile of this cluster indicates that the cluster 
may well spill out into two of the surrounding CCDs. To 
negate this problem we use only the furthest CCD from the 
cluster (with background area at least 12 arcmin distant) to 
estimate our field population. This solution, whilst increas- 
ing the statistical error on the blue fraction estimate, will 
prevent any cluster contamination of our field sample. 

MS2053. 7-0449: As previously noted in Section 4.2 we 
are not complete to a level M v = — 19 in this cluster. This 
is primarily due to the fact that some of our cluster images 
were contaminated with multiple tracks which we believe 
are reflected light from space debris. Once eliminated by 
masks designed for each individual image, this dropped the 
detection level below our initial target. We therefore avoid 
estimating the blue fraction at M v = —19, however our es- 
timate of f b at M v = —20 in (R — I) should still be valid. 
There may additionally be a low redshift cluster impinging 
upon MS2053.7 at about the R30 radius. Further calibra- 
tion images of this cluster are being obtained and extended 
results will be presented later in this series of papers. 

RXJ2106. 8-0510: This cluster is dominated by sev- 
eral large, bright elliptical galaxies, which in turn serve to 
obscure a large amount of the fainter, generally smaller, clus- 
ter population. The primary affect of this is that the esti- 
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Figure 5. Plot of blue fractions versus redshift for the eight clus- 
ters presented in this paper. The ft, results presented here are cal- 
culated using galaxies with magnitudes brighter than M v = —20, 
within R3o- Open circles (with error bars) represent ft, values for 
red colours (V — R and R — I for low and high redshift clus- 
ters respectively); diagonal crosses (with error ellipses) represent 
fb values for blue colours (B — V and V — R for low and high 
redshifts). The solid line represents the fit line from Butcher & 
Ocmlcr (1984), which should be compared with the open circles. 
RXJ2146.0 shows signs of foreground contamination (see text) 
but we include the results from this cluster for completeness, here 
represented by the solid circle (R — I) and solid square (V — R). 
We have also not included the V — R fraction from MS2053.7 as 
we are incomplete in the V band for this cluster (see text). The 
dot-outlined square represents the blue fraction range from the 
10 clusters in the X-ray selected sample of Smail et al. (1998). 
Error bars are one sigma, based on Poisson statistics only. 

mation of R30 is more troublesome than in other clusters. 
Additionally the CM diagram is less populated than in other 
cases, resulting in larger errors on the blue fraction. 

RXJ1606.7+2329 and RXJ1633.6+5714: Both 
these clusters are low X-ray luminosity poorer systems. This 
tends to increase the magnitude of the statistical error on the 
estimate of f b in both cases. Additionally the actual CMR 
is more difficult to fit (see Section 5.2). 
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Figure 6. Plot of blue fractions versus X-ray luminosity (0.5- 
2 keV) for the eight clusters presented in this paper. The f b re- 
sults presented here are calculated using galaxies with magnitudes 
brighter than M v = —20, within R3o- Symbols are the same as in 
Figure 5: open circles and vertical crosses (with error bars) rep- 
resent fb values measured in red colours; diagonal crosses (with 
error ellipses) represent ft, values measured in blue colours. Again 
RXJ2146.0 is represented by the solid symbols. 



6 DISCUSSION 

6.1 Blue fractions against redshift 

Wc find that most of our clusters have significant blue frac- 
tions. In Figure ^ we plot blue fraction versus cluster red- 
shift. The open circles represent V — R and R— I fb values for 
the low and high redshift clusters respectively; the diagonal 
crosses represent B — V and V — R fb values for the low and 
high redshift clusters. The solid square and circle represent 
the B — V and V — R values of RXJ2146. 0+0423 at z=0.532 
which suffers from contamination from a foreground group 
(see previous section), and will thus have excess blue galax- 
ies. We firstly consider the blue fractions measured from the 
reddest colours, which correspond approximately to B — V 
at rest (the circles in Figure |E|). Although a constant blue 
fraction of f b ~ 0.2 is a reasonable description of the current 
data, the blue fractions are also consistent with the trend- 
line of Butcher & Oemler (1984), which is overlayed on the 
data. A large scatter is also observed, although a large frac- 
tion of the scatter could in principle be due to measurement 
errors. The dot-outlined square represents the redshift and 
blue fraction (estimated using B-I colours) distribution of 
the 10 clusters in the Smail et al. (1998) sample. 

It is also interesting to note that in general there is not 
a good agreement between the results from the two colours. 
In nearly all cases the bluer colour gives higher fb values 
{ft ~ 0.4) than the red colour. The reasons for the ex- 
cess may be numerous. Firstly the intrinsic range in galaxy 
colours is larger in the bluest bands, as illustrated by the 
reduction in tightness of the CMR. This also has the effect 
of making the blue fraction cut criterion, estimated by the 
interpolation of predicted galaxy colours, more difficult to 
evaluate. Though we estimate that the colour of the blue 
galaxy cut would need to be lowered by at least 0.2 magni- 
tudes to decrease f b from 0.4 to 0.2 as measured in the red 



colour. This is significantly greater than any error on the 
interpolation between predicted galaxy colours. In fainter 
galaxies the scatter in galaxy colours, allied to the larger 
photometric errors, may contribute slightly to higher blue 
fractions. 

We believe, however, that this blue fraction excess can 
best be explained in terms of variations in the spectral en- 
ergy distributions (SEDs) of the galaxies. The larger range 
in the blue galaxy colours makes excess blueness in brighter 
galaxies (where the photometric errors are not dominant) 
much more apparent. The galaxies contributing to the ex- 
cess blue fraction often have B — V colours typical of Sab 
types whilst having colours indicative of later type galaxies 
in U — B (in general the galaxies that are most blue in one 
colour are also blue in the other). This is important when 
we consider that a fraction of the blue galaxies may have 
starburst or post-starburst SEDs (e.g. have the "E+A" or 
"k+a" spectra of many blue galaxies at high redshift (Couch 
& Sharpies 1987; Dressier et al. 1994; Couch et al. 1994)). In 
fact the observed B—R colours of the blue galaxies in our low 
redshift sample agree well with the Bj — Rf colours of the 
E+A galaxies found in Couch & Sharpies (1987), when filter 
conversion is accounted for. Spectroscopic observations will 
be needed to reliably determine the nature of these galaxies. 

One of the significant implications of the variation in 
blue fraction, as estimated from different colours, concerns 
sample selection. We argue that optical selection in any one 
filter may bias the cluster sample towards systems more ac- 
tive in the respective filter and conclude that a method in- 
dependent of the galaxies, such as X-ray selection, should 
be the preferred method of defining a truly representative 
sample. 

6.2 Blue fractions versus X-ray luminosity 

In Figure ^| we investigate the f b values as a function of X- 
ray luminosity. We conclude that no real trend can be found, 
over the factor of 100 or so in X-ray luminosity, though the 
large errors on the data points makes an expansion of the 
sample a priority before any rigorous conclusions can be 
drawn. This result is, however, in agreement with the lack of 
blue fraction correlation to X-ray luminosity in the Butcher 
& Oemler (1984) sample (Lea & Henry 1988; Andreon & 
Ettori 1999). If this result is substantiated, it has important 
implications for the physical processes that are responsible 
for the Butcher-Oemler effect. 

An interesting question is whether the galaxies within 
R30 would expect to have had their gas removed via ram- 
pressure stripping due to the ICM (Gunn & Gott 1972). 
An estimate of the dependence of ram-pressure stripping 
with L x can be made as follows. The rate of ram-pressure 
stripping is proportional to ~ Picm V ga i 2 ' 4 (Gaetz et al. 
1987) and L x oc Picm 2 ■ From the Virial theorem and the 
luminosity-temperature relation, V ga i oc T 0,5 oc Lx 1 / 6 , 
and so the ram-pressure stripping rate will increase as ~ 
L x . Thus the truncation of star-formation by ram-pressure 
stripping is expected to be significantly more effective in 
high L x clusters, especially over the factor of 100 in L x sam- 
pled here. 

Stevens et al. (1999) model a large spherical galaxy 
ploughing through a dense ICM and being stripped of gas. 
They find that gas stripping is primarily a function of clus- 
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Figure 7. Plot of blue fractions versus radius as a function of R30 for the eight clusters presented in this paper. The results shown here 
are for the V — R and R — I colours for the low and high redshift clusters respectively, corresponding approximately to B — V at rest. 
Clusters are shown in redshift order, from left to right, with low redshift clusters on the top row. The various symbols represent differing 
magnitude cuts. Open circles show blue fractions at M v < —21, diagonal crosses give ft, values for M v < —20 and triangles represent 
M v < —19. We give only results in the magnitude ranges in which we believe we are complete and at radii which are completely enclosed 
within the central CCD. Error bars are given at only one magnitude for clarity. X-ray luminosities are given in units of 10 44 ergs _1 . 
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Figure 8. Plot of blue fractions versus radius as a function of R30 for the eight clusters presented in this paper. The results shown 
here are for the B — V and V — R colours for the low and high redshift clusters respectively, corresponding approximately to U — B at 
rest. Symbols are as in Figure 7. 
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ter temperature (and hence mass), galaxy velocity within 
the ICM, the gas replenishment rate from stars and the dis- 
tance from the centre of the cluster potential. We can make 
an estimate of cluster mass from our X-ray luminosities 
via the X-ray luminosity-temperature relation (e.g. David 
et al. 1993) which has been shown to be non-evolving for the 
redshift ranges of our clusters (Mushotzky & Scharf 1997; 
Fairley et al. 2000). Our most massive cluster, MS1455.0, 
in which ram-pressure stripping should be most powerful, 
should have a temperature of approximately 6 keV. We esti- 
mate, from predictions presented in Figure 6 of Stevens et al. 
(1999), that ram-pressure stripping will only be a dominant 
force towards the centre of clusters, and then only for lumi- 
nous, high temperature clusters (i.e. T x > 3keV; see also 
Abadi et al. 1999). Thus, even if the orbits of most galaxies 
within R30 have taken them close to the cluster centre in 
the past, we would not expect ram-pressure stripping to be 
a dominant process in the low L x clusters. To increase the 
ram-pressure force, we would require that either infalling 
galaxies have velocities significantly larger than the mean 
galaxy velocity (which is not inconceivable, e.g. Henry & 
Lavery 1987), or that the gas replenishment rate is too low 
to be effective in replacing stripped gas, which may be less 
plausible. If the excess of blue galaxies in the low X-ray lu- 
minosity clusters studied here contains a large fraction of 
post-starburst galaxies, as found in richer clusters, then a 
mechanism other than ram-pressure stripping is probably 
responsible for truncating the star formation. 

6.3 Blue fractions with radius and galaxy infall 

It has been widely noted that blue fraction varies with 
cluster-centric radius. To investigate this in our sample we 
estimate the blue fractions detailed in the previous sections 
within various radii for each cluster. We consider here both 
colours and study the variation of fb with radius as a frac- 
tion of R30. Figures [7] and [] show plots of integrated blue 
fraction versus radius, for each of our clusters, in the red 
and blue colours respectively. We plot the blue fraction at 
varying magnitude cut-offs, in order to estimate the effect 
of including brighter or fainter galaxies. As can be clearly 
seen there is a large variation in the profiles of some of the 
clusters, a small change in others. In most of the clusters we 
see an increasing blue fraction with fainter magnitude cut 
(in agreement with Kodama & Bower (2001)). This trend, 
however, is not universal in our cluster sample. We note also 
that the slope of fb with radius, for each cluster, is similar in 
both colours, with the principle difference again being the 
magnitude of fb- 

Most of the clusters exhibit at least some increase in 
blue fraction with radius, which highlights the need for cau- 
tion when interpreting a blue fraction, at a characteristic 
radius, for a cluster. However the main conclusions drawn 
in the previous sections should not vary dramatically given 
small shifts in characteristic radii. The radius used here, R30, 
is of course dependent on the galaxy distribution. We also 
investigated the effect of using a radius derived from the X- 
ray properties of the cluster on the blue fractions. Here we 
used our estimated cluster temperatures to assign a cluster 
radius (r2oo), assuming the clusters followed an NFW profile 
(Navarro et al. 1995). We found that the blue fractions at 
r2oo/2 were in agreement (within errors) with those at R30. 



This is not too surprising as the galaxy population would be 
expected to roughly trace the cluster dark matter profile in 
relaxed clusters. 

To what extent are the conclusions drawn from Figures 
5 and 6 dependent on the absolute magnitude limit and the 
radius used to measure the blue fraction? Figures 7 & 8 
show that the broad trends of the variation of blue fraction 
with absolute magnitude, radius and colour are similar from 
cluster to cluster, so that the conclusions drawn from Figures 
5 and 6 are to a large degree independent of the precise 
choice of absolute magnitude limit or radius, so long as a 
consistent choice (with radius as a fraction of R30) is made. 

The increase in blue fraction with radius in many clus- 
ters suggests that it is field galaxies falling into the cluster 
which cause the enhanced blue population (e.g. Ellingson 
et al. 2000) and which may be one of the prime causes of the 
well known morphology-density relation, where spiral frac- 
tion increases with decreasing galaxy number density (e.g. 
Dressier et al. 1997). The mechanism by which the infalling 
field galaxies are morphologically transformed and have 
their star formation truncated could be ram-pressure strip- 
ping. As discussed in the previous section, the likelihood of 
ram-pressure stripping increases with ICM density towards 
the cluster centre. Gas stripping, however, should be less 
important in the lower luminosity clusters (e.g. RXJ1633.6 
and RXJ1606.7). Thus the increase in blue fraction with ra- 
dius seen in these systems should not be due to ram-pressure 
stripping. 



7 CONCLUSIONS 

We have presented the first results from a large multi-colour, 
wide-field imaging program of high redshift X-ray selected 
clusters. We use the large background area available from 
our observations of 8 clusters, using the INT WFC, to allow 
a statistical estimate of the fraction of blue galaxies in the 
cluster cores. The B, V, and R band photometry for four 
intermediate redshift clusters (z ~ 0.25) and V, R and I 
band photometry for four higher redshift systems (z ~ 0.5), 
corresponded approximately to UBV at rest. 

We find no significant change of blue fraction with red- 
shift in the range 0.2<z<0.5, but the blue fraction in all 
the clusters is higher than in those of Butcher & Oemler 
(1984) at z<0.1. Our reddest colour matches that used by 
Butcher & Oemler (1984). In this colour, the blue fractions 
we measure are equally consistent (within the measurement 
errors) with the trend-line presented in Butcher & Oem- 
ler (1984), or with no trend in blue fraction with redshift 
at 0.2<z<0.5. Our blue fractions in this colour are also in 
agreement with those measured in X-ray selected, luminous 
clusters by Kodama & Bower (2001) over a similar redshift 
range as used here. However, our blue fractions are higher 
than those measured in X-ray selected, luminous clusters 
by Smail et al. (1998) at zfa0.25. Smail et al. (1998) find a 
median blue fraction of fb = 0.04 ± 0.02 for concentrated 
clusters (C > 0.35). A possible explanation lies in the differ- 
ent colour (restframe U — Rot observed B — I) and selection 
band (restframe R) used by Smail et al. (1998) to define the 
blue fraction. 

We note a significant difference in the blue fraction val- 
ues as calculated in our two different colour bands. This may 
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explain some of the scatter in previous blue fraction mea- 
surements, which use a variety of different colours, and may 
additonally imply some level of bias in fb , depending on the 
optical band used for cluster selection. 

We find no evidence for a trend in blue fraction with 
X-ray luminosity, though again our results exhibit a large 
scatter. Interestingly we do not find the relation between 
blue fraction and richness that Margoniner et al. (2001) find. 
Unfortunately, whilst that study was based on a large clus- 
ter catalogue the authors chose to extract blue fractions at 
a common fixed spatial radius for each cluster. This obvi- 
ously prevents direct comparison with this work. It may also 
explain their richness correlation result, since a lower blue 
fraction in richer clusters could be due to the smaller radii 
(relative to R30) used to sample the richer clusters, combined 
with the morphology-density relation. 

Profiles of blue fractions with radius show an increasing 
proportion of blue galaxies, often faint, towards the outskirts 
of the clusters. The fact that the blue fraction is a function 
of radius, of faint magnitude cut, and also of the colour used 
to define it, emphasizes the need for extreme caution when 
assigning a characteristic blue fraction for an individual clus- 
ter. This may provide an explanation for the widely differ- 
ing blue fraction levels reported in various studies. Rakos & 
Schombert (1995), for example, find a high blue fraction of 
0.8 at z=0.9, but use a different method to most studies, 
employing Stromgren photometry, a variable cluster radius 
and photometric redshifts to help remove background galax- 
ies. Comparison of blue fractions measured using different 
techniques may not be valid. 

Our results suggest that the increased blue fractions 
are caused by infalling field galaxies, which may or may not 
undergo a starburst, before having their star formation trun- 
cated. We would expect that ram-pressure stripping occurs 
in the cores of our more X-ray luminous clusters, and this 
would cause the lowering of blue fraction, towards the clus- 
ter cores, that is observed. This explanation, however, is less 
likely to explain the same observational result in lower L x 
systems, where other physical processes (such as galaxy in- 
teractions) must be occurring. Our expanded cluster sample 
will allow a better understanding of the physical processes 
giving rise to blue galaxies in high redshift clusters. 
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